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Abstract 

Performing functional integration of the free Lagrangian, we find the vacuum energy of a field. 
The functional integration is performed in a way which easily generalizes to systems at non- 
zero temperature. We use this technique to obtain the Casimir energy density and pressure at 
arbitrary temperatures. 



1 Introduction 

If we have electromagnetic field in a restricted region, the vacuum energy of the electromagnetic field shows 
up as an effective zero point energy, which is known as Casimir energy. For a rectangular region bounded 
by two perfectly conducting plates separated by a distance £ along the z-axis, the magnitude of this energy 
is given by 

E = - W' (L1) 

where L x and L y are the dimensions of the plates (assumed large) along the other two axes. The energy 
density in the region between the plates is therefore given by 

755?- (L2) 

There are various ways that the magnitude of this energy density can be derived. Here, we outline a 
method using functional integration of the Lagrangian. We will show that this method easily extends to 
the case when the region between the two plates is at a finite temperature. Using this extension, we will 
calculate the Casimir energy at finite temperature. 

It is well known that the functional formulation of the electromagnetic field is complicated because of 
gauge invariance. The gauge volume has to be divided out by a Fadeev-Popov procedure. In order to avoid 
these complications, we perform all our calculation using a complex scalar field of mass M. In other words, 
we will calculate the zero-point energy of a complex scalar field in the specified geometry, assuming that 
this field satisfies the same boundary conditions as the electromagnetic field at the boundaries of the region. 
Once this is done, the result can be applied for the electromagnetic field as well, since it has the same number 
of degrees of freedom. We only need to put the mass of the field to be equal to zero. 



2 The functional formulation 



Let us then start from the Lagrangian of a free scalar field: 

C = (9V)(W-M 2 ^. (2.1) 
The generating functional obtained from this Lagrangian is given by 

Z = J [Zty][2tyt] exp (i J d A x£^j . (2.2) 

Since the free Lagrangian is quadratic, this integration can be performed formally, and one obtains 

Z = dct (d 2 + M 2 y 1 . (2.3) 

Thus, 

InZ = - tr In (<9 2 + M 2 ) , (2.4) 

using a well-known identity that for any operator, the logarithm of its determinant is the same as the trace 
of its logarithm. 

The trace mentioned above runs not only over any internal degrees of freedom that may be present (in 
the particular case at hand, there isn't any), but also over the space-time points. Thus, more explicitly we 
can write 

InZ = - J d*x (x\\n(d 2 + M 2 )\x) . (2.5) 
Going over to the momentum representation, this can be written as 

tlZ = ln H- 2 + M 2 ) 

where the integration over M 2 is indefinite. The integral on the right side is, of course, infinite. This is no 
surprise. For a general system, the quantity Z can be interpreted as 

lnZ = -z/ vac , (2.7) 

where / vac is the free energy density of the vacuum, so that we can write 

For a system at zero temperature, the free energy density is equal to the internal energy density, p vac - The 
infinity obtained for the integral on the right side then just corresponds to the well known result in quantum 
field theory that the vacuum energy of a field is formally infinite. Once normal ordered fields are employed, 
this infinite contribution vanishes. Effectively, this amounts to setting the zero of energy at the vacuum 
energy of the field. 

In the derivation above, we have implicitly assumed that the space-time region is infinite in all directions. 
For a more general situation, let us denote the obvious generalization of Eq. (2.8) by 

where the square brackets around the momentum integration measure now indicate that the measure has to 
be taken as appropriate in a particular situation. In fact, if the momentum values happen to be discrete, we 
should sum rather than integrate over the momenta. 

In order to calculate the vacuum energy at a non-zero temperature we can employ the imaginary 
time formalism, where the fco values for a bosonic field are quantized in the form 

k = 2win/@, (2.10) 
2 



for arbitrary integers n. In this case, the integration over ko is to be replaced by a sum over the imaginary 
energies. Thus, 



[d 4 fc] 

(2tt)4 - 
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[d 3 k] 



(2.11) 



We are still keeping the square brackets around the integration measure for the spatial components of 
momentum in view of the fact that in a restricted geometry as is relevant for the Casimir problem, the 
spatial components of momenta will also not be continuous. Thus, for a scalar field at finite temperature, 
we can write 
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[d 3 k] 



dM 2 
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13 n^ooJ ( 2?r ) 3 J (27rm//3) 2 - k 2 - M 2 

The sum over n present in Eq. ([2.12|) can be performed by using the identity 



oo 1 
X! 72 i 



i=-c 



I 2 +a 2 



■ coth 7ra . 



Writing 



coth z = 1 



a 2z 



we can rewrite Eq. (2.12) as 



[d 3 k] 

(2^)3 



dM 2 



2uj k 



using the shorthand 



L0 k 



. = Vk 2 + M 2 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



The integral over M 2 present in the resulting form is an indefinite integral. Instead of using k and M as 
independent parameters, we can use k and u>k ■ This change can be effected by making the substitution 



dM 2 = 2u>k duik ■ 



which gives 



fp 



[d 3 k] 
(2Vp 



dwh. 1 



( 1 + ^tt)-/- + ^ 



(2.17) 



(2.18) 



Clearly, the first term in the integrand gives a temperature independent part, which is the contribution at 
zero temperature (i.e., at (3 — > oo). As indicated, this term will be denoted by /oo. The second term is the 
temperature dependent contribution, which we have denoted by f'g. We now discuss these contributions one 
by one. 



3 The temperature independent contribution 

For the first term, which we denoted by , we can straightaway perform the indefinite integral over ujf. to 
write 



To proceed, we now use the identity 
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[d 3 k] 
J2^f 



(3.1) 



ds 

r(-i/2)7 

i r ds 
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exp [-s(k 2 + M 2 



(3.2) 



Substituting this form into the first term of Eq. (2.18), we can perform the fc integration if the space were 
infinite. This would give 



,-sM 2 



(3.3) 







ds 1 
s {Aits) 2 

This is essentially the form obtained by Schwinger || through his proper-time formalism. If the integration 
is performed now, we would get a result proportional to T(— 2), which would be infinite, as remarked earlier. 

But we also said that this energy is not really relevant for us. In the Casimir geometry, i.e., in the region 
between two infinite conducting plates at z = and z — £, the component of momentum perpendicular to 
the plates will be quantized. For an electromagnetic field, these quantized values will be given by 

k z = - (3.4) 

with arbitrary integers I so that the potential can vanish at both the plates. As we said in the Introduction, 
we will take this same boundary condition for the complex scalar field as well. Thus, for the Casimir 
geometry, we should write 
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(3.5) 



where k± indicates momentum values in the plane perpendicular to the z-axis. Using Eqs. (3T) and (3.2) 
as before, we can now perform the integration over k± to obtain 



/a 
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using the standard definition of the Jacobi 0-function[]: 
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(3.6) 



(3.7) 



The expression for f x given above can be simplified, noting that mathematically, the problem of Casimir 
geometry at zero temperature is equivalent to the problem of infinite geometry at non-zero temperature, 
tackled through the imaginary time formalism. We therefore follow a procedure enumerated by Dittrich f| 
for non-zero temperatures, adapting it to the present case. For this, one notes the following property of the 
Jacobi ^-function: 
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which enables us to write 
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Putting this form back in Eq. (3.6), we obtain 

ds 



(Atts) 
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(3.9) 



(3.10) 



There is a summation over I. Notice that, the I — term in this sum is exactly same as the result obtained 
in Eq. (^3) for infinite space. Thus, this result is also infinite. However, the physically important quantity 
is the difference of this quantity and the corresponding result for the infinite space, which is: 
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(3.11) 



See, e.g., Ref. 



(8.180, eq. 4 and the comment below. 



The integral over s can be expressed in terms of a modified Bessel function, defined byf] 

= * (2) / ^rexp . (3.12) 



Writing sM 2 = t in Eq. (|3.1lj) , the expression can be rewritten in the form 

2M 4 ™ [°° dt ( ( 2 l 2 M 2 \ M2 - 1 

This result is twice that of Ambjorn and Wolfram JjJ, who calculated the same quantity for a real scalar 
field. 

For the M — > limit which is relevant for finding the result for the electromagnetic field, we can use the 
form 

K 2 (z)^2z- 2 , (3.14) 

which is valid for small z. This gives 

This is the correct result for the boundary conditions explained above. As remarked in the Introduction, 
this result is applicable for the electromagnetic field which has the same number of degrees of freedom as 
the complex scalar field. Indeed, this result agrees with Eq. (1.2), since the free energy density in this case 
is equal to the internal energy density. 



4 The temperature dependence 

We now consider the temperature dependent part of the generating functional, which was given by the second 
term in Eq. ( [2.18 ) . Holstein and Pal M suggested that this term can be easily tackled if we rewrite it in the 



form 

00 



OC 



For infinite 3-dimensional space, the integral can directly be written in the form of a modified Bessel function 
H . To modify this procedure for the Casimir geometry, we first note that the exponent, written in terms of 
fc, has a square root owing to the definition of uik given in Eq. ( [2.16 ). This can be eliminated by using the 
identity 

lf°° ( z 2 \ 

e- = ^/ o ^^exp^--]. (4.2) 

Using this, for the Casimir geometry, we can utilize Eq. ( |3.5[ ) to write 

22., / r 2o2 2\ 
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Using 

u 2 = k\ + (ttI/C) 2 + M 2 (4.4) 

as is appropriate for this case, we can perform the integration over k± in a straight forward manner. This 
gives 

r* - -E, 2 U r***~*U-***\ ±e,J^l) 



V- 1 f; +1/2 ( r 2 (3 2 M 2 \ a f n nrr 2 [3 2 » 



2 See, e.g., Ref. [§, §8.432, eq. 6. 



Using the property of the i9-function given in Eq. (3.8), this can be rewritten in the form 



£ 



= -£ 



fri (^ 2 p 2 ) 2 j 

2 



dx x exp 



2 f} 2 M 2 
Ax 



Mo 



£ 



where, for the sake of notational simplicity, we have written 
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(4.7) 



The integral over x can be represented in terms of the modified Bessel function defined in Eq. fl3.12| ), and 
we obtain 



fa 
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J2 -js-K-iWMArt) . 



(4.8) 



For arbitrary values of M , this cannot be reduced further analytically. However, we arc interested in the 



limit M — » 0. From the general expression for modified Bessel functions given in Eq. ( 3.12 ), one can check 
that K u (z) = K_ v (z). Thus, we can write K 2 in place of K_ 2 in the last equation. Using now the limiting 



form of this function for small arguments which was given in Eq. (3.14), we obtain 
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h = ~ E 8^274 E 
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The sum over I can be performed by using the formula 
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which can be obtained by differentiating Eq. ( 2.13| ) with respect to the parameter a. Using this, we obtain 
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This result was obtained by various authors earlier, using different techniques §, |, fu| (y]|. 

One interesting check of this result is that, for £ — > oo, the expression in the square bracket has a limiting 
value of 2, so that 



#(c-oo) = -£- 



(ttt- 2 /? 2 ) 2 45/3 4 
The internal energy density can be calculated from here, which gives 

' f>+R 9f 'f> ^ 



df3 15/3 4 ' 



(4.12) 



(4.13) 



which is the well-known result for the energy density of a Planck distribution. 

Experimentally, however, we are interested about small values of Ci f° r which the remaining sum in 



Eq. (4.11) has to be performed numerically. The results are summarized in the next section. 



5 Discussion of the results 



Adding up the contributions given in Eqs. ( 3.15 ) and ( 4.11 ), we can write the total free energy density in 
the form 



720C 4 
G 



G(z), 



(5.1) 



where z is a dimensionless parameter defined as 



7T/3 



and 
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(5.2) 



(5.3) 



The corresponding density of the total internal energy will be given by 

7T 2 

720C 



08 = - ^KTlK(z) , (5.4) 



where 

K(z) = G(z)-zG'(z), (5.5) 
the prime on the function G implying differentiation with respect to its argument. This gives 



K{z) = 1 



E° z 4 ( r r / r \ 2 9 r / r \ 3 9 r r 1 , 
— ^-coth- + (-J csch 2 -+(-J csch 2 - coth - \ . (5.6) 

r—l K J 



Experimentally what is measured is the force on the condenser plates which define the boundaries of the 
region in the z-direction. The force per unit area, or the pressure p, is related to the total free energy F of 
a system through the thermodynamic relation 

—(5?).— 

where V is the volume of the region, and (3, or the temperature, is to be kept constant while taking the 
partial derivative. In the present situation, we can write this as 

fdfA 



\d( J 

Substituting the expression for fp obtained above, we obtain 



/' - f; -^^-\- (5 ' 8) 



240C 

where z is the variable defined in Eq. (|5.2|), and 



P=-^7I ff W' (5-9) 



H(z) = G{z)- l -zG'{z) 

OO 

= 1 - 30 V - csch 2 - coth - (5.10) 
z — ' r z z 

r=l 

Notice again that, in the limit C — > oo, this gives the familiar formula for the radiation pressure, which is 
one-third the internal energy density given in Eq. (4.13). 

The expressions for internal energy and pressure given above are valid if the region inside the capacitor 
plates are at a temperature 1/(3 and the outside region is at zero temperature. In an experiment, this is a 
difficult situation to realize, specially since the distance C between the plates is very small. Rather, a more 
practical scenario is when the region between the plates, as well as the ambient region, are both kept at a 
temperature 1/(3. In this case, the quantities of physical interest are the differences in the internal energy 
and pressure between the two regions — the one inside the plates, and the one outside. We denote these 
quantities by essentially the same notation used in Eqs. ( |5.4[ ) and (5.9, with two new functions K(z) and 
H(z) replacing K{z) and H (z) respectively. These new functions are given by 

K(z) = l-90V4!-coth-+ (-Y csch 2 - + (-Y csch 2 - coth -- 3) , (5.11) 
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Figure 1: The solid line is a plot of the function K (z), and the dotted line is of H(z). These functions give the 
difference between the internal energy density and pressure of the inside and outside regions. Both are normalized 
to have the value unity at zero temperature. 



and 



ff(z) = l-30]T 
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(5.12) 



In Fig. |l| we show these functions K(z) and H(z) as functions of the dimensionless variable z. The left end 
of the plots correspond to small z, i.e., large (3 or vanishing temperature. In this case, the results should 
reduce to the ordinary Casimir results, so that the gains should equal unity, as seen in the figure. As z 
increases, we see that the energy density quickly approaches the result appropriate for infinite volume. 

As far as the difference between the external and internal pressures is concerned, we see that the relevant 
function H(z) increases roughly linearly with z when z is large. Now, z can increase either because £ 
increases, or because the temperature increases. In the case when ( grows, the pressure difference, despite 



the increase of the function H(z), decreases owing to the factor of l/( 4 in Eq. (5.9), and vanishes in the 
infinite volume limit as expected. On the other hand, if we keep £ fixed and increase the temperature, we 
sec that although the difference in energy density between the external and the internal regions goes to zero, 
the pressure difference in fact increases with the temperature. This increase in the pressure difference is 
actually offset by an increase in entropy difference to keep the energy difference zero. 
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